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ABSTRACT. A parametric weighted graph is a graph whose edges are labeled with continuous real 
functions of a single common variable. For any instantiation of the variable, one obtains a standard 
edge-weighted graph. Parametric weighted graph problems are generalizations of weighted graph 
problems, and arise in various natural scenarios. Parametric weighted graph algorithms consist of 
two phases. A preprocessing phase whose input is a parametric weighted graph, and whose output 
is a data structure, the advice, that is later used by the instantiation phase, where a specific value for 
the variable is given. The instantiation phase outputs the solution to the (standard) weighted graph 
problem that arises from the instantiation. The goal is to have the running time of the instantiation 
phase supersede the running time of any algorithm that solves the weighted graph problem from 
scratch, by taking advantage of the advice. 

In this paper we construct several parametric algorithms for the shortest path problem. For the case 
of linear function weights we present an algorithm for the single source shortest path problem. Its 
preprocessing phase runs in 0(V 4 ) time, while its instantiation phase runs in only 0(E + V log V) 
time. The fastest standard algorithm for single source shortest path runs in 0(VE) time. For the case 
of weight functions defined by degree d polynomials, we present an algorithm with quasi-polynomial 
preprocessing time 0(V^ 1+log log v ) and instantiation time only 0(V). In fact, for any pair 
of vertices u, v, the instantiation phase computes the distance from u to v in only 0(log 2 V) time. 
Finally, for linear function weights, we present a randomized algorithm whose preprocessing time is 
0(V 3 ' 5 ) and so that for any pair of vertices u, v and any instantiation variable, the instantiation phase 
computes, in O(l) time, a length of a path from u to v that is at most (additively) e larger than the 
length of a shortest path. In particular, an all-pairs shortest path solution, up to an additive constant 
error, can be computed in 0(V 2 ) time. 
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1. Introduction 

In networking or telecommunications the search for the minimum-delay path (that is the short- 
est path between two points) is always on. The cost on each edge, that is the time taken for a signal 
to travel between two adjacent nodes of the network, is often a function of real time. Hence the 
shortest path between any two nodes changes with time. Of course one can run a shortest path al- 
gorithm every time a signal has to be sent, but usually some prior knowledge of the network graph 
is given in advance, such as the structure of the network graph and the cost functions on each edge 
(with time as a variable). 

How can one benefit from this extra information? One plausible way is to preprocess the 
initial information and store the preprocessed information. Every time the rest of the input is given, 
using the preprocessed information, one can solve the optimization problem faster than solving the 
problem from scratch. Even if the preprocessing step is expensive one would benefit by saving 
precious time whenever the optimal solution has to be computed. Also, if the same preprocessed 
information is used multiple times then the total amount of resources used will be less in the long 
run. 

Similar phenomena can be observed in various other combinatorial optimization problems that 
arise in practice; that is, a part of the input does not change with time and is known in advance. 
However, many times it is hard to make use of this extra information. 

In this paper we consider only those problems where the whole input is a weighted graph. 
We assume that the graph structure and some knowledge of how the weights on the edges are 
generated are known in advance. We call this the function-weighted graph - it is a graph whose 
edges are labeled with continuous real functions. When all the functions are univariate (and all 
have the same variable), the graph is called a parametric weighted graph. In other words, the graph 
is G = (V, E, W) where W : E — > F and F is the space of all real continuous functions with 
the variable x. If G is a parametric weighted graph, and r e lis any real number, then G(r) 
is the standard weighted graph where the weight of an edge e is defined to be (W(e))(r). We 
say that G(r) is an instantiation of G, since the variable x in each function is instantiated by the 
value r. Parametric weighted graphs are therefore, a generic instance of infinitely many instances 
of weighted graphs. 

The idea is to use the generic instance G to precompute some general generic information 1(G), 
such that for any given instantiation G(r), we will be able to use the precomputed information 1(G) 
in order to speed up the time to solve the given problem on G(r), faster than just solving the problem 
on G(r) from scratch. Let us make this notion more precise. 

A parametric weighted graph algorithm (or, for brevity, a parametric algorithm) consists of 
two phases. A preprocessing phase whose input is a parametric weighted graph G, and whose 
output is a data structure (the advice) that is later used by the instantiation phase, where a specific 
value r for the variable is given. The instantiation phase outputs the solution to the (standard) 
weighted graph problem on the weighted graph G(r). Naturally, the goal is to have the running 
time of the instantiation phase significantly smaller than the running time of any algorithm that 
solves the weighted graph problem from scratch, by taking advantage of the advice constructed in 
the preprocessing phase. Parametric algorithms are therefore evaluated by a pair of running times, 
the preprocessing time and the instantiation time. 

In this paper we show that parametric algorithms are beneficial for one of the most natural 
combinatorial optimization problems: the shortest path problem in directed graphs. Recall that 
given a directed real-weighted graph G, and two vertices u, v of G, the distance from u to v, denoted 
by 8(u, v), is the length of a shortest path from u to v. The single pair shortest path problem seeks to 
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compute S(u, v) and construct a shortest path from u to v. Likewise, the single source shortest path 
problem seeks to compute the distances and shortest paths from a given vertex to all other vertices, 
and the all pairs version seeks to compute distances and shortest paths between all ordered pairs of 
vertices. In some of our algorithms we forgo the calculation of the path itself to achieve a shorter 
instantiation time. In all those cases the algorithms can be easily modified to also output a shortest 
path, in which case their instantiation time is the sum of the time it takes to calculate the distance 
and a time linear in the size of the path to be output. 

Our first algorithm is a parametric algorithm for single source shortest path, in the case where 
the weights are linear functions. That is, each edge e is labeled with a function a e x + b e where a e 
and b e are reals. Such linear parametrization has practical importance. Indeed, in many problems 
the cost of an edge is composed from some constant term plus a term which is a factor of some 
commodity, whose cost varies (e.g. bank commissions, taxi fares, vehicle maintenance costs, and so 
on). Our parametric algorithm has preprocessing time 0(n 4 ) and instantiation time 0(m + n log n) 
(throughout this paper n and m denote the number of vertices and edges of a graph, respectively). 
We note that the fastest algorithm for the single source shortest path in real weighted directed graphs 
requires 0(nm) time; the Bellman-Ford algorithm Q. The idea of our preprocessing stage is to 
precompute some other linear functions, on the vertices, so that for every instantiation r, one can 
quickly determine whether G(r) has a negative cycle and otherwise use these functions to quickly 
produce a reweighing of the graph so as to obtain only nonnegative weights similar to the weights 
obtained by Johnson's algorithm lfl2l . In other words, we avoid the need to run the Bellman-Ford 
algorithm in the instantiation phase. The 0(n 4 ) time in the preprocessing phase comes from the use 
of Megiddo's[ 13] technique that we need in order to compute the linear vertex functions. 

Theorem 1.1. There exists a parametric algorithm for single source shortest path in graphs weighted 
by linear functions, whose preprocessing time is 0(n 4 ) and whose instantiation time is 0(m + 
n log n). 

Our next algorithm applies to a more general setting where the weights are polynomials of 
degree at most d. Furthermore, in this case our goal is to have the instantiation phase answering 
distance queries between any two vertices in sublinear time. Notice first that if we allow exponential 
preprocessing time, this goal can be easily achieved. This is not hard to see since the overall possible 
number of shortest paths (when x varies over the reals) is 0(n!), or from Fredman's decision tree 
for shortest paths whose height is 0(n 2 5 ) ©. But can we settle for sub-exponential preprocessing 
time and still be able to have sublinear instantiation time? Our next result achieves this goal. 

Theorem 1.2. There exists a parametric algorithm for the single pair shortest path problem in 
graphs weighted by degree d polynomials, whose preprocessing time is 0(n^°^ + ^° s ^ d ^ lo s n ) and 
instantiation time 0(log 2 n), where f(d) is the time required to compute the intersection points 
of two degree d polynomials. The size of the advice that the preprocessing algorithm produces is 

(9( n (0(l)+l°g d) logny 

The above result falls in the subject of sensitivity analysis where one is interested in studying 
the effect on the optimal solution as the value of the parameter changes. We give a linear-time 
(linear in the output size) algorithm that computes the breaking points. 

The practical and theoretical importance of shortest path problems lead several researchers to 
consider fast algorithms that settle for an approximate shortest path. For the general case (of real 
weighted digraphs) most of the algorithms guarantee an a-stretch factor. Namely, they compute a 
path whose length is at most a5(u, v). We mention here the (1 + e)-stretch algorithm of Zwick for 
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the all-pairs shortest path problem, that runs in 0(n w ) time when the weights are non-negative reals 
lTT8l . Here to < 2.376 is the matrix multiplication exponent Q. 

Here we consider probabilistic additive-approximation algorithms, or surplus algorithms, that 
work for linear weights which may have positive and negative values (as long as there is no negative 
weight cycle). We say that a shortest path algorithm has an e-surplus if it computes paths whose 
lengths are at most 5(u, v) + e. We are unaware of any truly subcubic algorithm that guarantees an 
e-surplus approximation, and which outperforms the fastest general all-pairs shortest path algorithm 
0. 

In the linear-parametric setting, it is easy to obtain e-surplus parametric algorithms whose pre- 
processing time is 0(n 4 ) time, and whose instantiation time, for any ordered pair of queried vertices 
u, v is constant. It is assumed instantiations are taken from some interval I whose length is inde- 
pendent of n. Indeed, we can partition I into 0(n) subintervals /j, J2, • • • of size 0(l/n) each, 
and solve, in cubic time (say, using 171'). the exact all-pairs solution for any instantiation r that is 
an endpoint of two consecutive intervals. Then, given any r G Ij = (aj,bj), we simply look at the 
solution for bj and notice that we are (additively) off from the right answer only by O(l). Standard 
scaling arguments can make the surplus smaller than e. But do we really need to spend 0(n 4 ) time 
for preprocessing? In other words, can we invest (significantly) less than 0(n 4 ) time and still be 
able to answer instantiated distance queries in 0(1) time? The following result gives a positive 
answer to this question. 

Theorem 1.3. Let e > 0, let [a, f3] be any fixed interval and let 7 be a fixed constant. Suppose G 
is a linear-parametric graph that has no negative weight cycles in the interval [a, /3], and for which 
every edge weight a e + xb e satisfies \a e \ < 7. There is a parametric randomized algorithm for the 
e-surplus shortest path problem, whose preprocessing time is 0(n 3 5 ) and whose instantiation time 
is 0(1) for a single pair, and hence 0(n 2 )for all pairs. 

We note that this algorithm works in the restricted addition-comparison model. We also note 
that given an ordered pair u, v and r G [a, /?], the algorithm outputs, in O(l) time, a weight of an 
actual path from u to v in G(r), and points to a linked list representing that path. Naturally, if one 
wants to output the vertices of this path then the time for this is linear in the length of the path. 

The rest of this paper is organized as follows. The next subsection shortly surveys related 
research on parametric shortest path problems. In the three sections following it we prove Theorems 
ll.llll.2l and ll.3l Section [5] contains some concluding remarks and open problems. 

1.1. Related research 

Several researchers have considered parametric versions of combinatorial optimization prob- 
lems. In particular function-weighted graphs (under different names) have been extensively studied 
in the subject of sensitivity analysis (see iPTTl ) where they study the effect on the optimal solution 
as the parameter value changes. 

Murty lfl4l showed that for parametric linear programming problems the optimal solution 
can change exponentially many times (exponential in the number of variables). Subsequently, 
Carstensen [3 ] has shown that there are constructions for which the number of shortest path changes 
while x varies over the reals is n n ^ ogn \ In fact, in her example each linear function is of the form 
a e + xb e and both a e and b e are positive, and x varies in [0, 00]. Carstensen also proved that this is 
tight. In other words, for any linear-parametric graph the number of changes in the shortest paths 
is n°( logn ). A simpler proof was obtained by Nikolova et al. ifToll . that also supply an n°( logn ' 
time algorithm to compute the path breakpoints. Their method, however, does not apply to the case 
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where the functions are not linear, such as in the case of degree d polynomials. Gusfield [ 10] also 
gave a proof for the upper bound of the number of breakpoints in the linear function version of the 
parametric shortest path problem, in addition to studying a number of other parametric problems. 

Karp and Orlin lfl5Tl . and, later, Young, Tarjan, and Orlin iMTl considered a special case of 
the linear-parametric shortest path problem. In their case, each edge weight e is either some fixed 
constant b e or is of the form b e — x. It is not too difficult to prove that for any given vertex v, when x 
varies from — oo to the largest xq for which G{x$) has no negative weight cycle (possibly xq = oo), 
then there are at most 0(n 2 ) distinct shortest path trees from v to all other vertices. Namely, for 
each r G [—00, xq] one of the trees in this family is a solution for single-source shortest path in 
G(r). The results in |[T5llT7l cleverly and compactly compute all these trees, and the latter does it 
in 0(nm + n 2 log n) time. 

2. Proof of Theorem O 

The proof of Theorem 11.1 I follows from the following two lemmas. 

Lemma 2.1. Given a linear-weighted graph G = (V, E, W), there exist a, ft £ MU{— oo}U{+oo} 
such that G{r) has no negative cycles if and only if a < r < (3. Moreover a and (3 can be found in 
0(n 4 ) time. 

Lemma 2.2. Let G = {V, E, W) be a linear-weighted graph. Also let a, (3 6KU {—00} U {+00} 
be such that at least one of them is finite and for all a > r > (3 the graph G{r) has no negative 
cycle. Then for every vertex v G V there exists a linear function g^'^ such that if the new weight 
function W' is given by 

W ((«, v)) = W ((«, v)) + gt'® ~ 9 [ * A 

then the new linear-weighted graph G' = (V, E, W) has the property that for any real a < r < f3 

all the edges in G'(r) are non-negative. Moreover the functions g^^ for all v G V can be found in 
0(mn) time. 

So given a linear-weighted graph G, we first use Lemma l2Tl to compute a and j3. If at least one 
of a and j3 is finite then using Lemma I2l2l we compute the n linear functions g^^\ one for each 
v G V. If a = —00 and j3 = +00, then using Lemma I2l2l we compute the 2n linear functions g^'°^ 
and <?l° ,/3 ' . These linear functions will be the advice that the preprocessing algorithm produces. The 
above lemmas guarantee us that the advice can be computed in time 0{n A ), that is the preprocessing 
time is 0(n 4 ). 

Now when computing the single source shortest path problem from vertex v for the graph G(r) 
our algorithm proceeds as follows: 

(1) If r < a or r > (3 output "—00" as there exists a negative cycle (such instances are consid- 
ered invalid). 

(2) If a < r < j3 and at least one of a or (3 is finite then compute g u (r) for all u G V. Use 
these to re-weight the edges in the graph as in Johnson's algorithm lTT2Tl . If a = —00 and 
f3 = +00 then if r < compute gjf (r) for all u G V and if r > compute gu'^ (r) for 
all u G V. Notice that after the reweighing we have an instance of G'(r). 

(3) Use Dijkstra's algorithm [6] to solve the single source shortest path problem in G'(r). Di- 
jkstra's algorithm applies since G'(r) has no negative weight edges. The shortest paths tree 
returned by Dijkstra's algorithms applied to G'(r) is also the shortest paths tree in G(r). As 
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in Johnson's algorithm, we use the results d'(v, u) of G'(r) to deduce d(v, u) in G(r) since, 

by Lemma l2.2l dfa. u) = d'(v, u) — g v (r) + g u (r). 
The running time of the instantiation phase is dominated by the running time of Dijkstra's algorithm 
which is 0(774 + n log n) 0. 

2.1. Proof of Lemma gH 

Since the weight on the edges of the graph G are linear functions, we have that the weight of 
any directed cycle in the graph is also a linear function. Let C\ , C2, ■ ■ ■ , Ct be the set of all directed 
cycles in the graph. The linear weight function of a cycle C, will be denoted by wt(Cj). If wt(Cj) 
is not the constant function, then let ji be the real number for which the linear equation wt(Cj) 
evaluates to 0. 

Let a and (3 be defined as follows: 

a = max{7j | wt(Cj) has a positive slope} . 

i 

(3 = min {73 | wt(Cj) has a negative slope} . 

i 

Note that if wt(Cj) has a positive slope then 73 = min^ {wt(Ci)(x) > 0} . Thus for all x > ji 
the value of wt(Cj) evaluated at x is non-negative. So by definition for all x > a the value of the 
wt(Cj) is non-negative if the slope of wt(Cj) is positive, and for any x < a there exists a cycle Cj 
such that wt(Cj) has positive slope and wt(C 4 )(x) is negative. Similarly, for all x < j3 the value of 
the wt(Cj) is non-negative if the slope of wt(Cj) is negative and for any x > (3 there exists a cycle 
Ci such that wt(Cj) has negative slope and wt(Cj)(x) is negative. 

This proves the existence of a and (3. There are, however, two bad cases that we wish to 
exclude. Notice that if a > P this means that for any evaluation at x, the resulting graph has a 
negative weight cycle. The same holds if there is some cycle for which wt(Cj) is constant and 
negative. Let us now show how a and f3 can be efficiently computed whenever these bad cases 
do not hold. Indeed, a is the solution to the following Linear Program (LP), which has a feasible 
solution if and only if the bad cases do not hold. 

Minimize x under the constraints 

Vi, wt(Ci)(x) > 0. 

This is an LP on one variable, but the number of constraints can be exponential. Using Megiddo's lfT3l 
technique for finding the minimum ratio cycles we can solve the linear-program in 0(n 4 log n) 
steps. 

2.2. Proof of Lemma H2] 

Let a and (3 be the two numbers such that for all a < r < (3 the graph G(r) has no negative 
cycles and at least one of a and f3 is finite. 

First let us consider the case when both a and f3 are finite. Recall that, given any number r, 
Johnson's algorithm associates a weight function h r : V — > M such that, for any edge (u, v) £ E, 

W {UtV) (r) + h r (u)-h r (v)>0. 
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(Johnson's algorithm computes this weight function by running the Bellman-Ford algorithm over 
G(r)). Define the weight function g%*'^ as 

This is actually the equation of the line joining (a, h a (v)) and (/?, h"{v)) in R 2 . 
Now we need to prove that for every a <r < (3 and for every (u, v) G V, 

W {u , v) (r)+gt'®(r)-9 [ y ] (r)>0. 

Since a < r < (3, one can write r = (1 — 5)a + 5f3 where 1 > 8 > 0. Then for all v G V, 

gfaflfr) = {1-S)h a {v) + Sh^(v) . 

Since W^ u , v ) ( r ) i s a linear function we can write 

W M (r) = (1 - <y)W (tt)t ,)(a) + 5W M (J3) . 

So after re-weighting the weight of the edge (u, u ) is 

(1 - S)W M (a) + W (u ,,)(/3) + (1 - S)h a (u) + <5^(u) - (1 - 8)h a (v) - 5h?(v) . 

Now this is non-negative as by the definition of hP and h a we know that both W( u ^ (/3) + h 13 (u) — 
h,P(v) and W( U)V ){a) + h a (u) — h a (v) are non-negative. 

We now consider the case when one of a or (3 is not finite. We will prove it for the case 
where f3 = +oo. The case a = — oo follows similarly. Consider the simple weighted graph 
Goo = (V, E, Woo) where the weight function is defined as: if the weight of the edge e is 
W(e) = a e x + b e then Woo(e) = a e . 

We run the Johnson's algorithm on the graph G^. Let h°°(v) denote the weight that Johnson's 

algorithm associates with the vertex v. Then define the weight function g^'°°^ as 

g^°°\x) = h a (v) + (x- a)h°°(v) . 
We need to prove that for every a < r and for every (u, v) G V, 

W M {r)+g^\r)-g^\r) = W M (r)+h a (u) + (r-a)h™ ( u )-h a (v)-(r-a)h°° (v) > 0. 

Let r = a + 6 where 5 > 0. By the linearity of W we can write W( u ^(r) = W( u ^{a) + 5a.( u>v y 
where W^ u , v ){ r ) = a (u,v) r + fyu.u)- So the above inequality can be restated as 

W M {a) + 5a M + h a (u) + 5h°°{u) - h a (v) - 5h°°{v) > . 

This now follows from the fact that both W( u ^{a) + h a {u) — h a (v) and a( M1) ) + h°°{u) — h°°(v) 
are non-negative. 

Since the running time of the reweighing part of Johnson's algorithm takes 0{mn) time, the 
overall running time of computing the functions g\?'^ is 0(mn), as claimed. 
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3. Proof of Theorem [L2] 

In this section we construct a parametric algorithm that computes the distance 6(u, v) between 
a given pair of vertices. If one is interested in the actual path realizing this distance, then it can be 
found with some extra book-keeping that we omit in the proof. 

The processing algorithm will output the following advice: for any pair (u,v) G V x V the 
advice consists of a set of t + 2 increasing real numbers — oo = &o < &i < • • • < b t < bt+i = oo 
and an ordered set of degree- d polynomials po, p\ , . . . , p t , such that for all b\ < r < bi+% the weight 
of a shortest path in G(r) from u to v is Pi(r). Note that each pi corresponds to the weight of a path 
from u to v. Thus if we are interested in computing the exact path then we need to keep track of the 
path corresponding to each pi. 

Given r, the instantiation algorithm has to find the i such that bi < r < bi + \ and then output 
Pi{r). So the output algorithm runs in time 0(log t). To prove our result we need to show that for 
any (u, v) G V x V we can find the advice in time 0(f(d)n) logn . In particular this will prove that 
t = 0(dn) logn and hence the result will follow. 

Definition 3.1. A minBase is a sequence of increasing real numbers — oo = bo < b\ < • • • < bt < 

bt+i = oo and an ordered set of degree-ci polynomials po, p\, . . . , p t , such that for all bi < r < bi + \ 
and all j / i,Pi(r) < pj(r). 

We call the sequence of real numbers the breaks. We call each interval [bi, &i+i] the i-th interval 
of the minBase and the polynomial pi the i-th polynomial. The size of the minBase is t. 

The final advice that the preprocessing algorithm produces is a minBase for every pair (u, v) G 
V x V where the i-th polynomial has the property that Pi(r) is the distance from u to v in G(r) for 
each bi < r < 

Definition 3.2. A minBase i (u, v) is a minBase corresponding to the ordered pair u, v, where the 
i-th polynomial pi has the property that for r G [bi, Pi(r) is the length of a shortest path from 
u to v in G(r), that is taken among all paths that use at most 2 l edges. 

A minBase e (u, w, v) is a minBase corresponding to the ordered triple (u, w, v) where the i-th 
polynomial pi has the property that for each r G p>i(r) is the sum of the lengths of a 

shortest path from u to w in G(r), among all paths that use at most 2 e edges, and a shortest path 
from w to v in G(r), among all paths that use at most 2 edges. 

Note that in both of the above definitions some of the polynomials can be +oo or — oo. 

Definition 3.3. If B\ and Bi are two minBases (not necessarily of the same size), with polynomials 
p\ and p?, we say that another minBase with breaks b' k and polynomials p' k is mvn{B\ + B^) if the 
following holds. 

(1) For all k there exist i, j such that p' k = pj + p 2 -, and 

(2) For b' k <r< b' k+1 and for all i, j we have p' k (r) <p\{r) +p 2 j{r). 

Definition 3.4. If B\, B2, . . . , B s are s minBases (not necessarily of the same size), with polynomi- 
als p\ ,p 2 , ■ ■ ■ ,pf , another minBase with breaks b' k and polynomials p' k is min{Bi, B%, . . . ,B S } 
if the following holds. 

(1) For all k there exist q such that p' k = pf , and 

(2) For b' k < r < b' k+1 and for all 1 < q < s and all i q , we have p' k {r) < p q iq {r). 
Note that using the above definition we can write the following two equations: 

minBase e+1 (u, v ) = min I minBase £ (u, w, v ) > . (3.1) 
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minBase e (u,w,v) = vain (minBase (u,w) + vainBase e (w,v) 



(3.2) 



The following claim will prove the result. The proof of the claim is omitted due to lack of 
space. 

Claim 3.5. If B\ and Bi are two minBases of sizes t\ and t 2 respectively, then 

(a) min(£>i + B 2 ) can be computed from B\ and B2 in time 0{t\ + t 2 )- 

(b) min{i?i, B2] can be computed from B\ and B2 in time 0(f(d){t\ + £2)), where f(d) is 
the time required to compute the intersection points of two degree-d polynomials. The size 
ofmin{£i,B 2 }isO(d(ti +t 2 )). 

In order to compute mm{B\, . . . , B s } one recursively computes X = min{i?i, . . . , B s / 2 } and 
Y = min{S s / 2+1 , . . . , B s } and then takes min{X, 1"}. 

If there are no negative cycles, then the advice that the instantiation algorithm needs from 
the preprocessing algorithm consists of vainBase^ ogn ^ (u,v). To deal with negative cycles, both 
minBase^ logn ^ (u, v) and minBase^ logn ^ +1 (u, v) are produced, and the instantiation algorithm 
compares them, if they are not equal, then the correct output is —00. 

Also note that minBase°(u, v) is the trivial minBase where the breaks are —00 and +00 and 
the polynomial is weight W((u, v)) associated to the edge (u, v) if (u, v) £ E and +00 otherwise. 

If the size of vainBase l {u, v) is Sf, then by (I3.ll ). (13 -2b . and by Claim 1331 the time to compute 
minBase e+1 (u,v) is 0{f{d)) Xogn si and the size of minBase i+1 (u,v) is 0(d) logn S£. Thus one 
can compute the advice for u and v in time 



and the length of the advice string is o(„(o(i)+i°g<9 i°g"). 

4. Proof of Theorem O 

Given the linear- weighted graph G = (V, E, W), our preprocessing phase begins by verifying 
that for all r € [a, /3], G(r) has no negative weight cycles. From the proof of Lemma I2T21 we know 
that this holds if and only if both G(a) and G((3) have no negative weight cycles. This, in turn, can 
be verified in 0(mn) time using the Bellman-Ford algorithm. We may now assume that G(r) has 
no negative cycles for any r G [a, /?]. Moreover, since our preprocessing algorithm will solve a large 
set of shortest path problems, each of them on a specific instantiation of G, we will first compute the 
reweighing functions of Lemma I2T21 which will enable us to apply, in some cases, algorithms 

that assume nonnegative edge weights. Recall that by Lemma [2721 the functions gl"'^ for all v G V 
are computed in 0(mn) time. 

The advice constructed by the preprocessing phase is composed of two distinct parts, which we 
respectively call the crude-short advice and the refined-long advice. We now describe each of them. 

For each edge e G E, the weight is a linear function w e = a e + xb e . Set K = 8(/3 — 
a) max e \ a e \. Let iVo = \Ky/n lnn/e] and let N\ = \Kn / e] . We define Nq + 1 and Ni + 1 points 
in [a, 0\ and solve certain variants of shortest path problems instantiated in these points. 

Consider first the case of splitting [a, f3] into iVo intervals. Let po = (f3 — a) /Nq and consider 
the points a + ipo for i = 0, . . . , Nq. The crude-short part of the preprocessing algorithm solves 
A^o + 1 limited all-pairs shortest path problems in G(a + ipo) for % = 0, . . . , Nq. Set t = 4-y/n In n, 
and let di(u, v) denote the length of a shortest path from u to v in G(a + ipo) that is chosen among 
all paths containing at most t vertices (possibly di(u, v) = 00 if no such path exists). Notice that 



(0(f(d)) logn ) logn 




176 



S. CHAKRABORTY, E. FISCHER, O. LACHISH, AND R. YUSTER 



di(u, v) is not necessarily the distance from u to v in G(a + ipo), since the latter may require more 
than t vertices. It is straightforward to compute shortest paths limited to at most k vertices (for 
any 1 < k < n) in a real- weighted directed graph with n vertices in time 0(n 3 log k) time, by 
the repeated squaring technique. In fact, they can be computed in 0(n 3 ) time (saving the log A; 
factor) using the method from HI, pp. 204-206. This algorithm also constructs the predecessor data 
structure that represents the actual paths. It follows that for each ordered pair of vertices u, v and 
for each i = 0, . . . , iVo, we can compute di(u, v) and apath j>j(u, v) yielding di(u, v) in G(a + ipo) 
in 0(n 3 |iVo|) time which is 0(n 3 5 Inn) . We also maintain, at no additional cost, linear functions 
fi(u, v) which sum the linear functions of the edges of pi(u, v). Note also that if di(u, v) = oo then 
Pi(u, v) and fi(u, v) are undefined. 

Consider next the case of splitting [a, (3] into Ni intervals. Let p\ = (/3 — a)/N\ and consider 
the points a + ip\ for i = 0, . . . , N\. However, unlike the crude-short part, the refined-long part 
of the preprocessing algorithm cannot afford to solve an all-pairs shortest path algorithm for each 
G(a + ipi), as the overall running time will be too large. Instead, we randomly select a set H C V 
of (at most) \fn vertices. H is constructed by performing ^fn independent trials, where in each 
trial, one vertex of V is chosen to H uniformly at random (notice that since the same vertex can be 
selected to H more than once \H\ < ^/n). For each h G H and for each i = 0, . . . , Nx, we solve 
the single source shortest path problem in G(a + ip±) from h, and also (by reversing the edges) 
solve the single-destination shortest path toward h. Notice that by using the reweighing functions 
g%*'® we can solve all of these single source problems using Dijkstra's algorithm. So, for all h G H 
and i = 0, . . . , Ni the overall running time is 

0(\Ni\\H\(m + nlogn)) = 0(n L5 m + n 2 ' 5 log n) = 0(n 35 ) . 

We therefore obtain, for each h G H and for each % = 0, . . . , N\, a shortest path tree Ti(h), together 
with distances d* (h, v) from h to each other vertex v G V, which is the distance from h to v in 
G(a + ipi). We also maintain the functions f*(h,v) that sum the linear equations on the path 
in T*(h) from h to v. Likewise, we obtain a "reversed" shortest path tree S*(h), together with 
distances d* (v, h) from each v G V to h, which is the distance from v to h in G(a + ipi). Similarly, 
we maintain the functions f*(v , h) that sum the linear equations on the path in S*(h) from v to h. 

Finally, for each ordered pair of vertices u, v and for each i = 0, . . . , N± we compute a vertex 
h u ,v,i G H which attains min^g^ d*(u, h) + d*(h, u) . Notice that the time to construct the h UjVj i 
for all ordered pairs u, v and for alH = 0, . . . , Ni is 0(n 3 5 ). This concludes the description of the 
preprocessing algorithm. Its overall runtime is thus 0(n 3 5 Inn). 

We now describe the instantiation phase. Given u, v G V and r G [a, f3] we proceed as follows. 
Let i be the index for which the number of the form a + ipo is closest to r. As we have the advice 
fi(u, v), we let wq = fi(u, v)(r) (recall that fi(u, v) is a function). Likewise, let j be the index for 
which the number of the form a + jpi is closest to r. As we have the advice h = h U)V j, we let 
^i = fj(u,h)(r) + f*(h,u)(r). Finally, our answer is z = min{wo, w\}. Clearly, the instantiation 
time is 0(1). Notice that if we also wish to output a path of weight z in Gir) we can easily do so 
by using either pi (u, v) , in the case where z = wq or using S* (h) and T* (h) (we take the path from 
u to h in S*(h) and concatenate it with the path from h to v in T*(h)) in the case where z = w\. 

It remains to show that, with very high probability, the result z that we obtain from the instanti- 
ation phase is at most e larger than the distance from u to v in G{r). For this purpose, we first need 
to prove that the random set H possesses some "hitting set" properties, with very high probability. 

For every pair of vertices u and v and parameter r, let p UjVjT be a shortest path in G(r) among 
all simple paths from u to v containing at least t = A^/n In n vertices (if G is strongly connected 
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then such a path always exist, and otherwise we can just put +00 for all u, v pairs for which no such 
path exists). The following simple lemma is used in an argument similar to one used in [ 18]. 

Lemma 4.1. For fixed u, v and r, with probability at least 1 — o(l/n 3 ) the path p u , v , r contains a 
vertex from H. 

Proof. Indeed, the path from p u , v ,r by its definition has at least Ay/nhin vertices. The probability 
that all of the ^fn independent selections to H failed to choose a vertex from this path is therefore 
at most 

< e" 41nn < — = o(l n 3 . 

n J n 4 

■ 

Let us return to the proof of Theorem 1 1.31 Suppose that the distance from u to v in G(r) is 5. 
We will prove that with probability 1 — o(l), H is such that for every u, v and r we have z < 5 + e 
(clearly z > 5 as it is the precise length of some path in G(r) from u to v). Assume first that 
there is a path p of length 5 in G(r) that uses less than 4-y/n In n edges. Consider the length of p in 
G[a + ipo). When going from r to a + ipo, each edge e with weight a e x + b e changed its length by 
at most |a e |po- By the definition of K, this is at most poK/(8((3 — a)). Thus, p changed its weight 
by at most 

It follows that the length of p in G(a + ipo) is less than 5 + e/2. But Pi(u, f ) is a shortest path from 
u to v in G(a + ipo) of ai l the paths that contain at most t vertices. In particular, di(u, v) < 5 + e/2. 
Consider the length of pi(u, v) in G(r). The same argument shows that the length of Pi(u, v) in 
G(r) changed by at most e/2. But wq = fi(u,v)(r) is that weight, and hence w < 5 + e. In 
particular, z < 5 + e. 

Assume next that every path of length 5 in G(r) uses at least In n edges. Let p be one such 
path. When going from r to r' = a + jp%, each edge e with weight a e x + b e changed its length by 
at most |a e |pi. By the definition of K, this is at most p\Kj (8(/3 — a)). Thus, p changed its weight 
by at most 

K K e 

In particular, the length of p UjV y is not more than the length of p in G(r'), which, in turn, is at 
most 5 + e/8. By Lemma |4~T1 with probability 1 — o(l/n 3 ), some vertex of h appears on p u ,v,r'- 
Moreover, by the union bound, with probability 1 — o(l) all paths of the type p u , v , r i (remember that 
r' can hold one of 0(n) possible values) are thus covered by the set H. Let h! be a vertex of H 
appearing in We therefore have d*Au, h') + d*(h', v) < 5 + e/8. Since h = h UjV ,j is taken 

as the vertex which minimizes these sums, we have, in particular, d*(u, h) + d*Ah, v) < S + e/8. 
Consider the path q in G(a + jpi) realizing d*(u, h) + d*(h, v). The same argument shows that 
the length of q in G(r) changed by at most e/8. But w\ = fj(u, h)(r) + fj{h, v)(r) is that weight, 
and hence w\ < 5 + e/4. In particular, z < 5 + e/4. 



5. Concluding remarks 

We have constructed several parametric shortest path algorithms, whose common feature is that 
they preprocess the generic instance and produce an advice that enables particular instantiations to 
be solved faster than running the standard weighted distance algorithm from scratch. It would be 
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of interest to improve upon any of these algorithms, either in their preprocessing time or in their 
instantiation time, or both. 

Perhaps the most challenging open problem is to improve the preprocessing time of Theorem 
ll.2l to a polynomial one, or, alternatively, prove an hardness result for this task. Perhaps less ambi- 
tious is the preprocessing time in Theorem ll.il 

Finally, parametric algorithms are of practical importance for other combinatorial optimization 
problems as well. It would be interesting to find applications where, indeed, a parametric algorithm 
can be truly beneficial, as it is in the case of shortest path problems. 
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